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Abstract 

The production of ir° and r\ mesons in collisions of light, intermediate-mass, and 
heavy nuclei at beam energies between 0.8A GeV and 2. OA GeV is discussed in the 
framework of a model assuming equilibrium at hadrochemical and thermal freeze- 
out, respectively. From the relative particle yields temperatures Tq of 56 - 90 MeV 
and baryon chemical potentials /ie of 780 - 675 MeV are deduced at chemical 
freeze-out. Midrapidity spectra of ir° and r\ mesons are consistent with the obtained 
temperatures if, in addition, the measured transverse expansion of the collision zone 
is taken into account. In contrast to results from similar analyses of AGS and SPS 
data the deduced freeze-out parameters are far below the phase boundary between 
a hadron gas and a quark-gluon plasma. 



1 Introduction 



During the last two decades the investigation of the properties of hadronic 
matter at high density and high temperature has become one of the key is- 
sues of nuclear physics. Relativistic nucleus-nucleus collisions offer the unique 
possibility to study nuclear matter under such extreme conditions. At inci- 
dent energies around 1A GeV the formation of a hot and dense reaction zone, 
the so-called nuclear fireball, has been verified experimentally [1]. According 
to model calculations, nuclear matter can be compressed for time spans of 
~10 - 15 fm/c to densities of ~2 - 3 times the nuclear groundstate density 
[2,3,4,5]. Simultaneously, temperatures up to ^100 MeV may be reached and 
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a substantial fraction of the nucleons participating in the collision is excited 
to heavier, short-lived resonance states which decay predominantly by meson 
emission [6,7,8]. Hence, the fireball at SIS-energies represents a system of in- 
teracting nucleons, resonances, and mesons. Under these conditions the onset 
of the partial restoration of chiral symmetry, which is a fundamental symmetry 
of quantum chromodynamics, could manifest itself in modifications of hadron 
properties inside the fireball. At beam energies of psIOA GeV, as they are avail- 
able at the AGS accelerator facility at BNL (Brookhaven), even higher baryon 
densities and temperatures are reached in the initial stage of the fireball [9]. In 
addition to the chiral phase transition one may approach in those collisions al- 
ready the phase transition from a hadron gas to a quark-gluon plasma due to 
the deconfinement of quarks and gluons. In such a scenario the fireball would 
no longer be purely hadronic during the whole collision process but the initial 
stage should be a partonic phase until hadronization takes place. At the SPS 
accelerator facility at CERN (Geneva), where the highest energy nuclear colli- 
sions can be studied using 158A GeV Pb beams, the possibility exists that the 
phase boundary to quark matter has been crossed already, although evidence 
for the phase transition is still rather tentative [10]. 

It is the common goal of numerous experiments at SIS, AGS, and SPS to 
investigate the hadronic or partonic fireball, respectively. Since each relativis- 
ts nucleus-nucleus collision, irrespective of the initial beam energy, develops 
a hadronic phase in its final stage, the importance of studying the purely 
hadronic fireball as it is formed in the SIS-energy regime is evident. Various 
observables, e.g. particle yields, spectra, and correlations are sensitive to dif- 
ferent parameters of the fireball. Only the combination of all observables may 
finally give a complete picture of A GeV heavy-ion collisions which up to now 
are not fully understood. 

In this paper we concentrate on ir° and 77-meson production as probes of 
the fireball. Chapter 2 summarizes briefly some basic features of the time 
evolution of relativistic nucleus-nucleus collisions as far as meson production 
is concerned. On the basis of the presently available systematics of 7r° and 
77-meson production the motivation for a thermal-model ansatz to describe 
hadronic matter at freeze-out is given. In the chapters 3 and 4 the model of a 
hadron gas in thermal and hadrochemical equilibrium is described in detail and 
applied to the fireball as observed at SIS. After fixing the model parameters 
by experimental data, the part of the nuclear-matter phase diagram probed 
by hadronic observables at SIS is determined and compared to corresponding 
results from AGS and SPS. Furthermore, the chemical composition of the 
fireball at freeze-out is estimated and the consistency of particle yields and 
spectra is studied. 
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2 7T° and r] Mesons from the Fireball at SIS 

2.1 Nucleus-Nucleus Collisions at SIS 



Numerous investigations of relativistic nucleus-nucleus reactions have shown 
that for central collisions at SIS the relative motion of target and projectile 
matter in the reaction zone almost comes to rest. While for heavy target- 
projectile combinations at moderate energies an isotropically expanding source 
at midrapidity has been observed [11], indications for incomplete stopping have 
been found in rapidity-density distributions for energies above ~ 1A GeV [12]. 
The size of the fireball thus arising depends on the centrality of the collision. 
Following the participant-spectator model, the geometrical overlap of the two 
nuclei determines the number of nucleons A part which are directly involved. In 
addition the energy available in the nucleon-nucleon system ^/s — 2m^ can 
be obtained from the projectile energy, so that the initial conditions of the 
system are fixed. In the subsequent development the energy is transformed 
into the excitation of thermal and collective degrees of freedom: it turns into 
heat and provides the mass stored in the resonance states, and it builds up 
compression and produces the flow of the expanding matter. 

In a simplified picture the mesons are trapped during the compression phase 
in a cyclic process of generation, absorption, and re-emission, exemplified for 
nucleons N, ir mesons, and A resonances by NN # NA # NNir. Within 
this approach the bulk of the mesons are released only with the onset of 
the expansion phase, when due to the decreasing matter density mesons and 
baryons decouple. This implies that hadronic probes reflect only the final state 
of nucleus-nucleus collisions and do not provide direct insight into the early, 
high-density phase of the collision as leptonic probes would do since they are 
not affected by the strong interaction. The moment in the expansion process 
when particles cease to interact strongly is called freeze-out. One distinguishes 
between chemical freeze-out and thermal freeze-out which correspond to those 
points in time when the relative abundances or the momentum distributions 
of the particles, respectively, are no longer changed. Obviously only inelastic 
collisions of nucleons, resonances, and mesons can change the relative particle 
yields while the momentum distributions are governed by the total interaction 
cross sections of these particles in the fireball. Therefore, chemical freeze- 
out does not occur after thermal freeze-out. Even in this simplified sudden 
freeze-out scenario one has to make sure that for the determination of freeze- 
out parameters only such observables are combined with each other which 
are fixed approximately simultaneously. In our ansatz particle yields reflect 
chemical freeze-out while momentum distributions reflect thermal freeze-out. 
The latter, however, are not purely thermal spectra in the sense that they 
follow Boltzmann distributions but they are modified by resonance decays 



3 



and reflect in addition collective flow phenomena. 

It is an interesting question to what extent the fireball at freeze-out can be 
described in terms of thermal and hadrochemical equilibrium. Equilibrium 
could be achieved through frequent scattering of particles in the fireball and it 
would be reached without any doubt if only the lifetime of the fireball would 
be sufficiently long. It has been shown that at AGS and SPS the hadronic ob- 
servables including strange particles are in quite good agreement with such a 
scenario [13,14,15,16], while at SIS the situation is still unclear. Only recently 
it has been reported that midrapidity momentum spectra of nonstrange par- 
ticles measured in central 58 Ni + 58 Ni collisions are consistent with thermal 
equilibrium in the energy regime from 1A GeV to 2A GeV if in addition col- 
lective radial flow is taken into account [12] . For the same reactions also agree- 
ment with chemical equilibrium has been claimed [17]. This analysis, however, 
was based on proton, deuteron, thermal charged pion, and A(1232)-resonance 
yields. The latter were obtained by decomposing the pion momentum spectra 
into a thermal component and into a component from resonance decay. In 
this paper we follow a different strategy: in order to keep the possibility to 
distinguish chemical and thermal freeze-out as individual stages in the time 
evolution of relativistic nucleus-nucleus collisions we strictly confine ourselves 
to an isolated discussion of relative particle yields and momentum distribu- 
tions, respectively, without combining these informations. In order to keep 
the thermal model as simple as possible we do not take into account parti- 
cles with open strangeness since at SIS energies their production rate is very 
small. Furthermore, an additional model parameter, the strangeness chem- 
ical potential /is, which, however, could be fixed by requesting strangeness 
conservation, would be necessary to incorporate strangeness consistently into 
the thermal model. With respect to strange particles, only recently Cleymans 
and co-workers have claimed that at 1.93 A GeV the particle yields including 
the strange hadrons agree with a complete chemical equilibrium scenario at 
freeze-out [18]. This would mean that not only the nonstrange and strange 
hadrons among themselves could be equilibrated seperately, but even equilib- 
rium between both the strange and the nonstrange sectors would have been 
established. 



2.2 Results from it and r] meson production experiments 

The neutral 7r° and r\ mesons which have large decay branches into pairs of pho- 
tons represent unique probes for the expanding fireball and therefore provide 
the basis for our discussion. Pions and T] mesons are the only mesons with sub- 
stantial production rates at SIS energies. While the pionic degree of freedom is 
also covered by the spectroscopy of the charged members of the multiplet, it is 
only through gamma-ray spectroscopy that the rj meson becomes observable 
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Fig. 1. Average 7r° and 77-meson multiplicities per average number of partici- 
pants in nucleus-nucleus collisions as a function of the energy available in the nu- 
cleon-nucleon system. The data are taken from [19,20,21,22,23,24,25,26,27,28,29]. 
The line represents a fit, taking into account also charged-meson production in 
addition [7]. 



in nucleus-nucleus collisions. These two neutral mesons can be simultaneously 
identified by a two-photon invariant-mass analysis of coincident photon pairs. 
During the last almost ten years the TAPS collaboration performed a series 
of systematic meson-production experiments covering incident energies from 
0.1 A GeV to 2.0 A GeV. These measurements have established an extended 
data base for n° and i] production in light ( 12 C + 12 C), intermediate ( 40 Ar, 

40 Ca + natQ^ 58 M + 58^ 86 Kr + nat^ and heayy systems (19? Au + 197 Au ) 

[19,20,21,22,23,24,25,26,27,28,29]. The primary information are the spectral 
distributions of the mesons, measured down to zero transverse momentum in 
narrow bins around midrapidity, and the meson yields, which are determined 
in the same rapidity intervals and are then extrapolated to the full solid angle. 
All data concerning the experiments 40 Ar + nat C& at 0.18 A GeV [27], 58 Ni + 
58 Ni at 1.9 A GeV [24], and 40 Ca + na *Ca at 2.0 A GeV [28] presented in this 
paper are not yet submitted for publication and should therefore be regarded 
as preliminary results. 

The first important result is obtained from the systematic investigation of the 
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meson yields. Fig. 1 shows the average meson multiplicities < M >, measured 
in inclusive reactions of various systems, normalized to the average number 
of participating nucleons < A part > as a function of the energy ^/s — 2itin 
available in the NN system normalized to the corresponding meson mass. A 
steep rise of the relative yields with increasing energy is visible. In addition, 
a weaker dependence on the size of the collision system is observed with the 
clear tendency for smaller inclusive yields in the heavier systems. This has 
also been observed in charged-pion production experiments [30,31]. In first 
approximation all data points fall onto a smooth curve indicating that the 
meson-production probabilities mainly depend on the available energy. This 
is quite remarkable because 7r° and rj production proceed through different 
baryon resonances: ir° mesons mainly originate from A(1232)-resonance de- 
cays, while the heavier 77 mesons essentially come from the N(1535) resonance 
which, at SIS energies, is the only significantly populated baryon resonance 
with a large decay channel into 77 mesons. The fact that the specific production 
mechanism seems not to influence the observed meson yields significantly may 
be interpreted as a first indication for meson emission from an equilibrated 
source. 

The second important result concerns the average transverse momenta < pt >. 
Fig. 2 shows < p t > of tt° and 77 mesons as a function of the energy available 
in the NN system. Although the experimental uncertainties for the 77 data 
are large compared to the 7T° data, the 77 transverse-momentum spectra reveal 
at a given energy a tendency towards larger mean values < p t > which is 
manifest around 2. OA GeV beam energy. To give a measure for the expected 
order of magnitude of < p t >, it can be noted that for the decay of A(1232) 
and iV(1535) resonances at rest one would get < p t > 7T o= 225 MeVQ and 
< Pt >r)= 182 MeV, respectively. For both mesons < p t > increases with 
the available energy which could be interpreted in terms of a considerable 
heating of the collision zone or a gradually growing population of the full 
resonance-mass distributions. For tt° mesons the rise of the average transverse 
momenta becomes smaller with increasing energy which is supported by a 
measurement of < p t > for tt~ mesons produced at 3.37^4 GeV bombarding 
energy [32]. The steep rise of the meson-production yields on the one hand 
and the onset of a saturation of average transverse momenta on the other 
hand clearly indicate that with increasing bombarding energy the fraction of 
available energy converted into massive particles becomes larger compared to 
the chaotic motion of particles in the fireball. 

Finally, the distributions of the transverse mass m t of tt° and 77 mesons are 
discussed. At midrapidity the transverse mass of a particle is equivalent to its 
total energy in the center-of-mass system. Following [33] the transverse-mass 
distribution of particles with mass m, emitted isotropically from a thermal 



Throughout this paper we follow the usual convention c = h = 1. 
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Fig. 2. Average ir and 77-meson transverse momenta measured in narrow rapid- 
ity intervals around midrapidity as a function of the energy available in the nu- 
cleon-nucleon system. The data are taken from the same experiments as listed in 
Fig. 1. The data points around available energies of 365 MeV (0.8^4 GeV beam en- 
ergy) and 447 MeV (LOA GeV beam energy) are shifted slightly on the energy axis 
with respect to their nominal position to resolve the error bars. Lines interpolating 
the 7T° and 77 data, respectively, are drawn to guide the eyes. A ir~ point taken from 
[32] has been added to the data set. 

source and characterized by a Boltzmann temperature Tg, can approximately 
be described by 

1 da (— m t \ . , / " ^ 

— - — oc exp — — with m t = Jm 2 + pf (1) 

mf dm t \ T B J 



In Fig. 3 transverse-mass distributions of 7r° and 77 mesons are plotted together 
with fits according to the parametrization given in Eq. 1 for three different 
collision systems at incident energies near 2A GeV. Within each reaction the 
7T° and 77 spectra exhibit almost identical inverse-slope parameters Tb which in 
addition do not change significantly with the mass of the colliding nuclei. Such 
a behaviour, the so-called m t scaling, is well known in high-energy physics 
[34,35]. Having in mind the different production mechanisms for both meson 
species it is furthermore quite remarkable that for m t > the 7r° and r\ 
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Fig. 3. Transverse-mass spectra of ir° and r/-mesons as observed in the systems 
12 C+ 12 C [25] and 40 Ca+ na *Ca [28] at 2.0A GeV beam energy and in 58 Ni+ 58 Ni 
[24] at 1.9A GeV. The distributions are divided by the square of the transverse 
mass. In this representation midrapidity particle spectra from a source in thermal 
equilibrium are expected to exhibit a purely exponential spectrum. The solid lines 
represent exponential fits to the 7r° data for m t > 400 MeV. 

data roughly fall on one commen line for all systems. This indicates that the 
energy required to produce a given transverse mass completely determines 
the relative abundance of the meson species near midrapidity. For low m t , 
however, individual differences become apparent. Fig. 3 shows a systematic 
enhancement over the exponential rise extrapolated from the high-m t region 
if one goes from the light to the heavier systems. Possible mechanisms that 
have been suggested to account for this fact involve enhanced pion rescattering 
through resonance states in the heavier systems [36] and multiple pion decay 
of heavy resonances [37]. The same observation, namely the m t scaling on 
an absolute scale and the low-mi enhancement in heavier systems, is not 
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only manifest at 2A GeV but has been recognized at various energies down to 
0.8A GeV and seems to be a general feature of heavy-ion collisions in the SIS— 
energy regime [20,24,25,26,28,29]. Although m t scaling does not prove that the 
meson source actually has reached thermal or hadrochemical equilibrium, the 
particle spectra of such an equilibrated source would be dominated by phase 
space as it is observed. 

In summary, the systematic study of 7T° and r) meson yields and momentum 
distributions reveals strong indications for the hypothesis that these mesons 
may be emitted from a source in hadrochemical and thermal equilibrium. 



3 Particle Yields and Chemical Freeze— Out in a Hadron— Gas Model 

3.1 Formalism of the Hadron-Gas Model 

In our ansatz we assume that at the point of chemical freeze-out the fireball 
can be described in terms of an equilibrated hadron gas. As constituents of 
this ideal gas we take into account nucleons, deuterons, pions, r\ mesons, and 
all nonstrange baryon resonances up to a mass of 1.76 GeV, which corresponds 
to yfs — ni>N for the nucleon-nucleon system at 2 A GeV beam energy. 

Nuclear matter in hadrochemical equilibrium has a particle composition char- 
acterized by a temperature T c and a baryon chemical potential Consider- 
ing the grand-canonical ensemble description of noninteracting fermions and 
bosons the particle-number densities pi are given by integrals over the particle 
momentum p: 



where gi is the spin-isospin degeneracy, Ei the total energy in the local rest- 
frame, and Bi the baryon number of the particle species %. The ± sign refers 
to the different statistics for fermions or bosons, respectively. 

Eq. 2 can not be applied to the fireball at SIS directly but only after the in- 
clusion of two modifications. To account for the fact that the fireball occupies 
only a finite volume, a surface correction has to be included. For an estimate 
of this correction we assume a spherical volume with radius R giving a mo- 
mentum dependent correction factor f(pR) [38,39] which has to be multiplied 
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with the integrand in Eq. 2: 



f(pR) = 1 



3tt 1 



(3) 



ApR ( P R) 2 



Introducing this correction leads to a decrease of individual particle-number 
densities down to 60 - 70 % of the uncorrected values for reasonable radius 
parameters. Particle yield ratios, however, are much less affected by the cor- 
rection since there the correction almost cancels out. The radius parameter 
should be chosen such that the number of baryons inside the spherical volume 
is in reasonable agreement with the average number of nucleons participating 
in the collision. We use a constant radius parameter of 5.5 fm for all beam en- 
ergies. If one assumes that chemical freeze-out may occur at baryon densities 
in the range between 0.2 and 0.6 po (see below) the corresponding freeze-out 
volume would then contain ^23 - 70 baryons. We verified that in our model 
analysis a variation of R by ±2 fm changes the resulting temperatures and 
chemical potentials by less than 5 %. 

The second modification is related to the fact that resonances do not have a 
fixed mass but exhibit a broad mass distribution. In contrast to thermal-model 
analyses of AGS and SPS data, a reasonable treatment of mass distributions 
can not be neglected in the 1A GeV energy regime since in this case the 
energy available in the NN system is small compared to the difference in mass 
between the nucleon and the resonances except for the A(1232) resonance. 
That is the reason why the A(1232) is by far the most abundantly populated 
resonance and why all the other resonances are populated predominantly in 
the tails of their mass distributions towards lower masses. Eq. 2 therefore has 
to be folded with normalized mass distributions A^m) which are 5 functions 
for all particle species except resonances. The following scheme to parametrize 
resonance mass distributions is to a large extent adopted from [37]. 

We distribute the resonance masses according to Lorentzian functions A R {m) 
which are determined by the mean resonance masses m R and the mass depen- 
dent total and partial decay widths T(m): 



We use the following parametrizations for the partial decay widths T(m) where 
the mean masses m R and widths Y R at m R are taken from [40]: 

• For the pion decay width of the A(1232) resonance we adopt the parametriza- 
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tion given by Koch et al. [41] 



where m is the actual mass of the A(1232) resonance, q and q R are the 




pion three momenta in the restframe of the resonance with mass m and 
rriA, respectively. The parameter 5 in the cutoff function has the value 5 = 



• The l7r and rj decay widths of the higher baryon resonances are given by 



where Bi WjV is the branching ratio of the resonance into the lir or 77 decay 
channel at m R , respectively, / is the angular momentum of the emitted 
meson, and q and q R are the momenta of the meson in the restframe of the 
decaying resonance as defined above. In this case we use 



It has to be noted that the only resonance decaying into the r\ channel in 
our ansatz is the iV(1535) resonance. 

• The 2-7T decay widths of higher baryon resonances are described in terms of 
one-step processes. This is in contrast to [37] where in a first step a higher 
baryon resonance decays into a A(1232) or iV(1440) resonance and a pion 
or into a nucleon and a p or a meson, and where in a second step the in- 
termediate resonance decays after propagation through the nuclear medium 
into a nucleon and a pion or into a nucleon and two pions. Since the popu- 
lation of higher resonances is small compared to the A(1232) population we 
simplify this procedure by treating the 2n decay as a direct process where 
the resonance decays into a nucleon and an object with angular momentum 
I = and twice the pion mass which subsequently decays into two pions. 
The partial decay widths can then be written in complete analogy to the 
In and r\ decay width. 

• Additional decay channels are neglected. 

Furthermore, one could consider an excluded-volume correction to take into 
account the hadron-hadron hard-core repulsion or in other words the differ- 
ence between a real and an ideal hadron gas. Since such a correction does not 
play a significant role for particle-yield ratios as demonstrated in [16] we stick 
to the ideal hadron-gas model. Another argument to omit the excluded vol- 
ume correction is the fact that chemical freeze-out seems to occur at baryon 
densities of the order of half the nuclear groundstate density or even less as 
discussed below. The hadron gas therefore is rather dilute which justifies the 
assumption of non-interacting particles at chemical freeze-out. 



300 MeV. 
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Finally, the isospin asymmetry in collisions of heavy nuclei, which has to be 
considered if the model analysis for example relies on charged-pion yields, 
is not taken into account in the present analysis. This is adequate since our 
approach is based on the isospin-independent r^-meson yields and the neutral- 
pion yields which in the isobar model correspond to one third of the total pion 
yields irrespective of the actual isospin asymmetry in the system. 

In summary, the equation we finally use to evaluate the particle-number den- 
sities within our hadron-gas model reads as follows: 



9* f tt t>\ *a f Ai(m)dm 
An J J exp \[ym 2 + p 2 - fi B ^i) /T c 



± 1 



(8) 



Because we have parametrized the resonance-mass distributions using mass- 
dependent partial decay widths we can easily calculate the particle-number 
densities of mesons originating from any individual resonance decay, which 
turns out to be essential for the determination of the model parameters from 
experimental data. If p[^ 2 n,r] * s defined as the particle-number density of pions 
or ^-mesons stemming from a specific decay channel of an individual resonance 
one obtains 



oo 



271 { { r tot (m)exp [(^m 2 +p 2 - ^ B B t ) T c \ ± 1 



with w — 1 in the case of one-pion or r^-meson decay and w = 2 for the 
two-pion decay of a resonance. 



3.2 Determination of the Model Parameters from Experimental Data 



For a given beam energy the values of the two parameters of the hadron-gas 
model - the temperature Tc and the baryon chemical potential hb - can be 
derived from at least two relative abundances of the constituents. The 7r° and 
?7-meson yields are particularly suited for such an analysis, as these particles 
unambiguously arise from the fireball. We base our analysis on the system- 
size averaged ratios < M n o > / < A part > and < M rj > / < M n o > as 
compiled in Tab. 1 (see also Fig. 1). The measured ratios are related to the 
particle-number densities calculated within the model by two equations: 




'p,+ /3 

v -Res / 

PN + 2PD + E PRes 
Res 



(Pb,T c ) 



(10) 
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Table 1 

Mean inclusive ir° multiplicities < M n o > relative to the average number of par- 
ticipating nucleons < A part > and mean inclusive jy-meson multiplicities < M v > 
relative to < M n o > measured for various target-projectile combinations in the 
beam-energy range from 0.8A GeV to 2.0^4 GeV. The system-size averaged values 
are given in addition. The 2.0A GeV Ca and 1.9A GeV Ni data are still preliminary. 
*For the average the latter have been upscaled by «10 % to 2. OA GeV using the fit 
of the energy dependence as shown in Fig. 1. 



E [A GeV] 


System 


<M jr o> [o/i 


<m; > [%] 


Reference 


0.8 


12 C+ 12 C 


3.7 ± 0.3 


0.31 ± 0.11 


[25] 


0.8 


40 Ar+ na *Ca 


3.1 ± 0.5 


0.41 ± 0.04 


[26] 


0.8 


197 Au+ 197 Au 


1.6 ± 0.3 


0.38 ± 0.10 


[29] 


0.8 


1 O < 1 (17 i 

12 C + 197 Au 


1.8 ± 0.2 




[25] 


0.8 


197 Au + 12 C 


2.0 ± 0.2 




[25] 


0.8 


Average 


2.6 ± 0.3 


0.39 ± 0.06 




1.0 


12 C+ 12 C 


5.6 ± 0.4 


0.57 ± 0.14 


[25] 


1.0 


40 Ar+ na *Ca 


3.0 ± 0.3 


1.3 ± 0.8 


[20,21] 


1.0 


^Kr+ raa *Zr 


2.5 ± 0.7 


1.3 ± 0.5 


[20,21,22] 


1.0 


197 Au+ 197 Au 


2.3 ± 0.5 


1.4 ± 0.5 


[20,21,22] 


1.0 


Average 


4.0 ± 0.4 


1.06 ± 0.41 




1.5 


40 Ar+ na *Ca 


6.5 ± 0.5 


2.2 ± 0.4 


[20,21] 


2.0 


12 C+ 12 C 


13.8 ± 1.4 


3.6 ± 0.4 


[25] 


2.0 


40 Ca+ na *Ca 


11.1 ± 1.1 


2.7 ± 0.3 


[28] 


1.9 


58 Ni+ 58 M 


8.6 ± 0.9 


3.3 ± 0.4 


[24] 


2.0 


Average* 


11.3 ± 1.1 


3.34 ± 0.37 





" ^ " 1 /;/ (HB,T C ) (11) 



exp 



Ptt 

Res 



Epfr + p% 8 ) /3 

Res / 



To deduce [Ib and T c we determine within our model the chemical potentials 
and temperatures for which the calculated ratios are in agreement with the 
experimental values for < M n a > / < A part > and < M v > / < M n a >, 
respectively. Each ratio therefore defines an area in the parameter plane Tc 
versus Fig. 4 shows the resulting areas for the four considered beam ener- 
gies 0.8, 1.0, 1.5, and 2.0A GeV.. The ratio < M^o > / < A part > is strongly 
correlated with fi B . To increase < M n o > / < A part > smaller values for 
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Fig. 4. Determination of chemical freeze-out parameters: for a given beam energy 
the two measured particle ratios < M n o > / < A part > and < M v > / < M n o > 
define two bands in the parameter plane (Tc/ie) as described in the text. The 
four panels show these bands for beam energies of 0.8, 1.0, 1.5, and 2.0A GeV. The 
freeze-out parameters are determined by the Tc and \i b values where the two bands 
cross each other. Except for the 0.8^4 GeV case additional bands are shown which 
are defined by deuteron-to-nucleon ratios measured in central 58 Ni+ 58 Ni collisions 
near the considered beam energies [12]. 

the chemical potential of the hadron gas are needed. In contrast to that, the 
ratio < M v > / < M n o > is mainly correlated with Tc- With increasing 
< M v > I < M w o > larger values for the hadron-gas temperature are nec- 
essary to obtain agreement between experimental and calculated ratios. For 
all beam energies the areas defined by the experimental boundary conditions 
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Table 2 

Chemical freeze-out parameters Tq and \x b obtained from inclusive 7r° and 77-meson 
yields for the four considered beam energies. 



E [A GeV] 


T c [MeV] 


fi B [MeV] 


0.8 


56 ± 2 


780 ± 10 


1.0 


70 ± 5 


770 ± 15 


1.5 


83 ± 7 


745 ± 20 


2.0 


90 ± 6 


675 ± 20 



cross each other exactly once and therefore determine the temperature and 
the baryon chemical potential of the hadron gas at chemical freeze-out in an 
unambiguous way. In addition, bands corresponding to deuteron-to-nucleon 
ratios measured by the FOPI collaboration [12] in central 58 Ni + 58 Ni colli- 
sions at beam energies of 1.06A GeV, 1A5A GeV, and 1.93A GeV, respectively, 
are shown in Fig. 4. Obviously also these ratios are in nice agreement with 
the chemical freeze-out parameters obtained from inclusive it and 77-meson 
yields. The results are summarized in Tab. 2. The analysis reveals a system- 
atic increase of the freeze-out temperature from 56 to 90 MeV between 0.8 
and 2.0A GeV which is accompanied by a reduction of the baryon chemical 
potential from 780 to 675 MeV. The corresponding baryon density at chemical 
freeze-out is approximately 0.2 - 0.6 times the nuclear groundstate density 
which is in reasonable agreement with results from other analyses even at 
higher energies [14,15,17]. 

Fig. 5 shows our results within the phase diagram of hadronic matter together 
with other (Tc, Hb) values obtained from particle-production experiments at 
SIS, AGS, and SPS energies [14,15,17,42]. The data points are connected by a 
line to symbolize the chemical freeze-out curve of hadronic matter. In contrast 
to the AGS and SPS results the deduced chemical freeze-out parameters at 
SIS energies are far below the phase boundary between a hadron gas and a 
quark-gluon plasma. Results dealing with thermal freeze-out, which are also 
shown in Fig. 5, are discussed in the following chapter. 

Having determined the chemical freeze-out parameters one can study the 
chemical composition of the hadronic fireball within the hadron-gas model. 
Fig. 6 shows the calculated relative populations of baryon species at chem- 
ical freeze-out as a function of the energy available in the nucleon-nucleon 
system together with the system-size averaged relative meson multiplicities 
< M n o yrj > I < A part > as given in Tab. I. The hatched areas correspond to 
the relative populations of nucleons and deuterons (summed), A(1232) and 
A(1535) resonances as the main sources of pions and rj mesons, respectively, 
and the sum of all remaining A and N resonances. With increasing available 
energy the fraction of baryons excited into resonance states at chemical freeze- 
out grows from ^2 % at 0.8A GeV beam energy to «15 % at 2.0A GeV. The 
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Fig. 5. Phase diagram of hadronic matter: the temperatures T and baryon chemi- 
cal potentials [is derived from yield ratios of particles produced in nucleus-nucleus 
collisions at different incident energies. Results of this work are plotted together 
with results taken from [14,15,17,42]. The solid curve through these data points 
represents the curve for chemical freeze-out of hadronic matter. Open data points 
represent parameter pairs for thermal freeze-out and are connected by a correspond- 
ing freeze-out curve (dashed). 

A(1232) resonance is populated most abundantly since it is the lowest mass 
baryon resonance. The ratio of heavier resonances to the A(1232) increases 
from s*5 % to ~15 % if one goes from low to high energies. Furthermore, it 
is obvious that the decay of resonances still excited at chemical freeze-out is 
not sufficient to account for the meson multiplicities observed asymptotically 
in the detector. Consequently, a sizeable fraction of resonance decays occurs 
before freeze-out is reached: at chemical freeze-out ~60 % of the pions and 
~90 % of the 7] mesons are not bound in baryon resonances but are already 
present as free mesons. 

Finally, one can extrapolate the resonance population at chemical freeze-out 
to the population in the high-density phase of the collision. Guided by mi- 
croscopic model calculations [37,43] which quite consistently give a maximum 
baryon density of p max ~ 2.5 p an d which predict ratios around 0.6 for the 
number of A(1232) resonances present at maximum compression relative to 
the number of asymptotically observed pions, one obtains an initial resonance 
density of 0.6 x 0.113 x 3 ~ 0.2 p max « 0.5 po at 2A GeV beam energy. Since in 
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Fig. 6. Inclusive 7r° and 77-meson yields as function of the energy available in the 
nucleon-nucleon system (left scale). The data points represent system-size averaged 
values from various collision systems (see Tab. 1). In addition, the hatched bands 
show the relative populations of nucleons and deuterons (summed), the A (1232) 
resonance, the iV(1535) resonance, and the sum of all remaining A and N resonances 
(right scale), as obtained from the hadron-gas model described in the text. Only 
one third of the A(1232) resonances decays into neutral pions. Considering that due 
to the limited available energy the iV(1535) resonance is predominantly populated 
at masses below the nominal mass of 1535 MeV it is comprehensible that less than 
50 %, which is the branching ratio at the nominal mass, of these resonances decay 
into the r] channel. 

the high-density phase of the collision the baryon-resonance density becomes 
comparable to the nuclear groundstate density, hadronic matter in that stage 
is often referred to as resonance matter [6,8]. 
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4 Particle Spectra and Thermal Freeze— Out in a Thermal Model 



While particle yields contain information about chemical freeze-out parame- 
ters as discussed in the preceeding chapter, particle spectra reflect the nucleus- 
nucleus collision process at the stage of thermal freeze-out in our simple sud- 
den freeze-out scenario. Even if the system reaches global thermal equilibrium 
at freeze-out the spectra, however, do not follow pure Boltzmann distribution 
functions, but they are modified by resonance decays and they furthermore 
reflect collective flow phenomena. 

We present here an analysis of the midrapidity 7r° and 77-meson spectra within 
the blast model proposed by Siemens and Rasmussen [44]. In this model, 
a fireball in thermal equilibrium expands isotropically. All particles in the 
fireball share a common temperature T and radial-flow velocity f3 r , which are 
the two parameters that characterize the fireball at thermal freeze-out. The 
transverse-mass spectra are then described by 



1 d 2 N 
mf dm t dy 



oc cosh?/exp(— •y r E/T) 



T\ sinha 



a 



— cosh a 
E 



(12) 



where 7 r = — (3%, a = (jrPrP) /T, and E, p, and y are the total en- 

ergy, momentum and rapidity of the considered particle in the center-of-mass 
system. The modification of the spectra compared to pure Boltzmann distri- 
butions becomes more significant the heavier the considered particle species 
is. 

Since pions and rj mesons have low masses their spectra are not very sensi- 
tive to the radial-flow velocity. Therefore a fit to the transverse-mass spectra 
using T and f3 r both as free parameters allows to determine the freeze-out 
parameters only with very large uncertainties. To overcome this difficulty, we 
take f3 r from the systematics of transverse-flow velocities for different collision 
systems and bombarding energies as published in [45]. Values for f3 r between 
0.25 and 0.35 are found to be a reasonable choice in the considered range of 
incident beam energies. To determine the thermal freeze-out parameter T we 
then fit the measured midrapidity transverse-mass spectra of 7r° and r\ mesons 
using Eq. 12, where we vary only T and keep f3 r constant. In the case of the 
pion spectra, we furthermore restrict the fit to m t > 400 MeV since the low m t 
region is expected to be modified by resonance decays. As an example Fig. 7 
shows the meson m t spectra measured in the systems 40 Ar/ 40 Ca + ncrf Ca to- 
gether with two fits for each considered beam energy. The fits corresponds to 
values of f3 r = 0.25 and f3 r = 0.35, respectively. Obviously, the distributions are 
consistent with a thermal freeze-out scenario as assumed in the blast model. 
This is also the case for the other investigated projectile-target combinations. 
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Fig. 7. Transverse-mass spectra of ir° and rj mesons measured in the system 
40 Ar/ 40 Ca+ na *Ca [20,21,26,28] at beam energies of 0.8, 1.0, 1.5, and 2.0,4 GeV. 
The cross sections are scaled in arbitrary units. For each beam energy two lines are 
drawn which represent fits to the 7r° data above transverse masses of 400 MeV using 
Eq. 12 with fixed transverse-expansion velocities of (3 r = 0.25 and 0.35, respectively. 
On the scale chosen a difference between the two fits is not visible. 

Fig. 8 shows the thermal freeze-out temperatures, again averaged over system 
size, as a function of the incident energy, together with the chemical freeze-out 
temperatures as obtained in the previous chapter. In addition, temperatures 
obtained in a blast-model analysis of 7r~, proton, and deuteron spectra from 
the FOPI collaboration [12] are shown for comparision. With increasing beam 
energy the thermal as well as the chemical freeze-out temperatures grow. 
Within the error bars, chemical and thermal freeze-out temperatures overlap 
and agree nicely with the FOPI results. Our analysis, however, shows some 
trend towards lower temperatures at thermal freeze-out which one could ex- 
pect since thermal freeze-out should not occur before chemical freeze-out. 
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Fig. 8. System-size averaged thermal and chemical freeze-out temperatures, respec- 
tively, as function of the beam energy. Thermal freeze-out temperatures taken from 
[12] are shown for comparision. 



Finally, we want to include our results on thermal freeze-out parameters into 
the phase diagram of hadronic matter as shown in Fig 5. To this aim one 
has to estimate the baryon chemical potential /ig at the moment of thermal 
freeze-out. If one assumes an isentropic expansion of the collision zone one 
can make use of the fact that the ratio T//j,b remains constant in that case. 
The baryon chemical potential at thermal freeze-out may then be calculated 
from the one at chemical freeze-out just by multiplication with the ratio of 
the temperatures. The resulting thermal freeze-out parameters are included in 
Fig. 5 together with a curve representing the thermal freeze-out of hadronic 
matter. While at SIS energies thermal and chemical freeze-out seem to al- 
most coincide, at higher energies the freeze-out parameters clearly indicate 
that thermal freeze-out reflects a later stage of the collision process as lower 
temperatures compared to chemical freeze-out are deduced. 
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5 Conclusion 



Neutral-pion and 77-meson data measured near midrapidity in the energy range 
from 0.8A GeV to 2.0A GeV are consistent with the assumption that the 
hadronic fireball, which evolves during relativistic nucleus-nucleus collisions, 
may be in hadrochemical or thermal equilibrium at the moment of chemical 
or thermal freeze-out, respectively. 

In contrast to ultrarelativistic collisions studied at AGS and SPS, chemical and 
thermal freeze-out seem to occur at almost the same stage of the expansion 
process at SIS energies. Another difference is that the freeze-out parameters 
at SIS are far below the phase boundary between a hadron gas and a quark- 
gluon plasma. This implies that the early stage of an ultrarelativistic nucleus- 
nucleus collision may come close to that phase transition while at SIS the 
early, high-density phase is described best by the term resonance matter. 

Neutral-meson yields and spectra reflect only one facet of the very complex 
process of a relativistic nucleus-nucleus collision. With these observable alone 
it is not possible to decide whether chemical or thermal equilibrium, respec- 
tively, is actually reached during any stage of such a collision. Further ex- 
periments and analyses investigating complementary probes are necessary to 
clarify to what extent nucleus-nucleus collisions in the SIS-energy range can 
be described within the framework of equilibrium models. 
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